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Abstract 

To date it has not been possible to prove whether or not the three-dimensional 
incompressible Euler equations develop singular behaviour in finite time. Some pos- 
sible singular scenarios, as for instance shock- waves, are very important from a 
physical point of view, since they imply the connection among the macroscopic 
and the microscopic scale. Therefore, the appearence of this type of singularity 
or a similar one might be interpreted as a possible explanation of the transition 
to turbulence. In order to clarify the question, some one-dimensional models for 
ideal incompressible hydrodynamics have been introduced and analyzed, and it was 
proven that shock-waves appear in finite time within this type of flow. In this work 
we question the validity of these models and analyze the physical meaning that the 
occurrence of a singularity in an incompressible flow, if it happens, may have. 
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1 Introduction 



It is not yet known whether or not the three-dimensional incompressible Eu- 
ler equations develop singular behaviour in finite time. Far from being a pure 
mathematical problem, it has deep connections with some fundamental ques- 
tions in physics. Actually, the relation between the appearance of finite time 
singularities and the onset of turbulence has been conjectured, because singu- 
larity formation may be a mechanism of energy transfer from large to small 
scales (or may not, see [1,2]). There are several types of blow-ups that can oc- 
curr in a fiuid fiow, such as, for instance, an infinite increment of the velocity 
in some spatial point. This kind of blow-up is not useful for an understanding 
of turbulence, since it might be interpreted as a break down of the nonrela- 
tivistic description. This is, if the velocity becomes too large, it is necessary 
to modify the equations to include relativistic effects [3] . On the other hand, 
a shock-wave is a singularity of the first spatial derivative of the velocity. 
Shock-waves are events in the fiow where the continuum description breaks 
down and a kinetic description, such as the Boltzmann equation, becomes 
necessary [4,5]. This suggests that shock- waves are a link between the micro- 
scopic and the macroscopic scale. In conclusion, we may argue that the onset 
of turbulence might be related with such situations in which the fiow develops 
a shock- wave (or a similar phenomenon, see below), but in which the velocity 
remains bounded. All along this work we will be concerned (unless explicitely 
indicated) with fiows carrying a finite total kinetic energy, since we consider 
this characteristic as necessary to keep the physical meaning of the fiow. 

First of all, let us clarify a question of language. Usually, it is understood that 
a shock-wave is a discontinuity of the velocity of the fiow. Of course, the first 
spatial derivative of the fiow will present a Dirac delta singularity located at 
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the shock, but it is possible to find a singularity of the first spatial derivative 
of a continuous velocity field, just by letting the slope of the tangent go to 
infinity. We can name this type of events (continuous velocity and divergence 
of its first spatial derivative) as quasi-shock-waves, due to its high similarity 
to traditional shock-waves. Let us explain why. Consider a continuous one- 
dimensional velocity field v{x), its first spatial derivative at xq is given by 
v'{xq) = limft^o('i^(^o + h) — v{xo))/h. This quantity will be divergent if there 
is a high accumulation of different values of v near xq. A high accumulation of 
different values of the velocity in the neighbourhood of one point reveals that 
the continuum description of the fluid is less accurate, and this descripcion 
becomes ill defined when a divergence is developed. This shows that shock- 
waves and quasi-shock-waves are singular events of the very same kind, both 
revealing the failure of the continuum description of the fiuid fiow and the need 
of a kinetic description (and maybe the connection between the macroscopic 
and the microscopic scales). An equivalent way to describe the phenomenon 
is to see it as the collision of two fiuid particles, as we will see below. 

2 First approach to singular behaviour 

The three-dimensional incompressible Euler equations read: 



where v is the velocity, a three-dimensional vector, and p is the pressure, a 
scalar. Eq.(la) is Newton's law for the fiuid, the left hand side is the convective 
derivative of the velocity and the right hand side is the force in terms of the 



dtv + {v ■ V)v = — Vp, 



(la) 



V • t; = 0, 



(lb) 
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pressure. Eq.(lb) is the incompressibility condition for the fluid. These equa- 
tions can be obtained from the Boltzmann equation performing the Chapman- 
Enskog expansion [6], but full mathematical proofs are available only before 
the appearance of singularities [7,8]. The Boltzmann equation describes the 
dynamics of a rarefled gas, taking into account two basic processes: the free 
flight of the particles and their collisions. In the hydrodynamic description, the 
convective derivative simulates the free flight of the particles and the gradient 
of the pressure is the flrst order correction to the interaction among them. The 
second order correction yields the Navier-Stokes equations, that are obtained 
by adding the term fiV'^v to the right hand side of Eq.(la), where /j, > 
is the viscosity of the fluid. The third order is still an open problem, how- 
ever some insight has been gained into it by means of a regularization of the 
Chapman- Enskog expansion [9,10]. These partial achivements suggest that at 
small scales the viscosity loses effectivity [10], what makes more interesting 
to study the dynamics of an inviscid fluid when we are close to the onset of 
turbulence. 

Due to the mathematical difficulties that a direct treatment of the Euler equa- 
tions imply, many simplifled models have been introduced in order to un- 
derstand better Euler dynamics. One classical approximation is the Burgers 
equation: 

dtv + vd^v = 0, (2) 

where v is the one-dimensional velocity. It is interesting to note that Burgers 
equation is a one-dimensional analogue of Euler equations where the pressure 
has been supressed. This equation is known to develop shock-waves in flnite 
time as a consequence of the crossing of its characteristics [11]. This result 
can be easily understood if one remembers the form of the Hamilton-Jacobi 



4 



equation for the free particle 




(3) 



recalling that the momentum is given hy p — d^S we see that deriving once 
with respect to the spatial variable Eq.(3), substituting p — mv, and reversing 
the direction of time t ^ —t (remember that classical mechanics is invariant 
with respect to time reversing) we recover Burgers equation. Thus Burgers 
equation can be thought as the evolution of the velocity of a cloud of free 
particles. Since this gas is noninteracting, the characteristics of the Hamilton- 
Jacobi equation, the Hamilton equations (the particle trajectories), will cross 
when subject to an adequate initial condition, what physically means that the 
fluid particles collide in finite time. Obviously, this fact does not require an 
infinite kinetic energy or an unbounded velocity of the fiuid. This cannot be 
directly translated to the Euler equations, since the pressure term acts as a 
force against this type of collisions. However, if we introduce an interaction 
between the particles in the form of a viscosity 



we get global existence of the solution. So it seems that in this case the shocks 
are an artifact of the noninteracting character of the particles. 

It is important to understand what we mean by fiuid particles and by collisions 
among them. A fiuid particle is not a molecule composing the fluid. It is a 
coarse-graining concept that corresponds to a point in the hydrodynamical 
(macroscopic) description of the fiuid. But it is actually composed of a high 
number of molecules. Correspondingly, a collision among fluid particles does 
not implie the collapse of two molecules. It is only the superposition of two 
coarse-graining entities reflecting the failure of the macroscopic approach. It 
is then necessary to describe the fluid microscopically in terms of the actual 



(4) 
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molecules composing the flow (this is the kinetic decription we refer to above) . 



It is known that the compressible Euler equation develop shock-waves as 
well [12] (and that they are again regularized by the viscosity in one dimen- 
sion). This shocks are, however, completely due to the compressible character 
of the flow, and cannot be extrapolated to incompressible flows [12]. This sug- 
gests that, if a singularity is present in an incompressible flow, it would not be 
of the shock type. The real problem is in contrast much more complex: when 
small-scale structures appear through the nonlinear dynamical evolution, they 
tend to display, at least locally, a much faster dependence on one particular 
spatial dimension (a phenomenon known as depletion of nonlinearity) [13]. 
In the limit, a flow collapsed into one spatial dimension displays a singular 
behaviour of the quasi-shock-wave type: the fluid particles collapse into a one- 
dimensional structure revealing the failure of the continuum description. Also, 
there is rigorous work showing that if a singularity is present, it is due to the 
collapse of a small-scale structure [14] . This shows that if a divergence appears 
in the three-dimensional Euler flow, it might be a phenomenon related to a 
shock-wave, although not necessarily a shock-wave itself. 

To clarify the question at hand, let us consider one example of infinite energy 
exact solution that shows a strong tendence towards the shock-wave. Some 
exact infinite energy singular solutions might be found in [15], and among 
them we choose the following one in cylindrical coordinates 



V = v{vr{r, t),0,v^{z,t)), 



(5) 



without a swirling component of the velocity, and where 



Vr{r, t) 



1 r 
2T* - 1' 



(6a) 
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v,{z,t) = - y, (6b) 

where T* is the blow-up time. For a better understanding of the fluid dynamics 

one can integrate the equations of motion for the fluid particles 

dr . , ,_ , 

— ^Vr{r,t), (7a) 



dz 

^ = ^.(^,^), (7b) 



that yield 



^it)=^i^Nj^f (8a) 



9{t)=9{0), (8b) 

m-m'^- (8c) 
This clearly indicates that all the particles of the fluid collapse in the plane 
z = ai the flnite time t = T*. In spite of this strong tendence to form a 
shock, this type of blow-up cannot be considered as a shock-wave since the 
collapse occurs in a spatial point at an inflnite distance from the origin, and 
at an infinite velocity. What this solution indicates is that the nonrelativistic 
description of the flow has broken down. 



3 Model equations for the vorticity dynamics 

One of the most important results proved about the regularity of the solutions 
of the Euler equations is the Beale-Kato-Majda theorem [16], that says that 
the solution exists globally in time if and only if 

/ max\uj( X, s)\ds < oo, (9) 

JO ^ 
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where u — V x v is the vorticity. This made very interesting to study the 
evolution of the vorticity in the Euler equations: 



dtuj + {v ■ V)uj ^u-Vv, (10) 
where the velocity can be recovered from the Biot-Savart Law: 

This equation can also be expressed in the following way 

V{uj)uj, (12) 



Dt 



where Duj/Dt = dtUJ + {v ■ V)a; is the convective derivative of the vorticity, 
and P is a symmetric matrix given by 



V = (A,) = 



1 / dvi_ ^ dvj_ 

2 \dxj dxi 



(13) 



The operator relating uj to Vuj is a linear singular integral operator that com- 
mutes with translation. In one spatial dimension there is only one such oper- 
ator: the Hilbert transform. 

In this spirit, Constantin et al. [17] proposed the following one-dimensional 
model for the vorticity equation: 

dtUJ = H{uj)uj, (14) 

where H{ijj) is the Hilbert transform of uj: 

H(.) = ip.V. r ^dy, (15) 

TT J-oo X — y 

and P.V. denotes the Cauchy principal value integral 

/oo r 
f{x)dx = lim J ^ f{x)dx. (16) 
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This equation was solved explicitly and it was shown that it blows up for some 
finite time Tq. However, it has been proven that the viscous analogue of this 
equation 

dtou = H{uj)uj + eujxx, (17) 
blows up for some finite time such that Tg < Tq; this is, adding diffusion 
makes the solution less regular. This is unsatisfactory in view of the result 
by Constantin [18], which says that if the solution to the Euler equations is 
smooth then the solutions to the slightly viscous Navier-Stokes equations are 
also smooth. In order to prevent this behaviour, De Gregorio [19,20] introduced 
an improved model keeping the convective derivative: 

dtUJ + udxU) = ujHiuj) + iidxx^ (18) 

with viscosity /i > 0. This equation does not develop singular behaviour, and 
De Gregorio concluded that one-dimensional models for hydrodynamics are 
not able to faithfully represent three-dimensional incompressible fiow. 

Baker, Li, and Morlet studied a very similar one-dimensional model simulating 
vortex sheet dynamics [21]: 

dtO - H{e)dj = lie,,, (19) 

which has been reinterpreted as a model for the Euler equations by Cordoba, 
Cordoba, and Fontelos [22] after switching = 0. In this model ^ is a scalar 
carried by the flow, and the vorticity and the velocity are defined, respectively, 
hy cu — 9x and v — —H{9). In reference [22] it is proven that this system blows 
up for some finite time in some spatial point, provided the initial condition 
is even, positive and compactly supported; then limj^T. ||6'a;||L°° = oo for 
some T* < oo. As proven there, the solution is even whenever the initial 
condition is so, implying that its Hilbert transform, say the velocity, is odd, 
and thus null at the origin. Also, the transport character of this equation 
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implies that ||^||Loc.(i) = ||^||loo(0). Numerical evidence has indicated that the 
blow-up appears in the origin [22,23], where a cusp of 6 is formed. This type 
of singularity corresponds to a (quasi-) shock, but we still need to know if the 
velocity is bounded at the origin to assure that a genuine (quasi-) shock- wave is 
formed. It is an important fact that the positivity of 9 implies that v — —H{6) 
is a decreasing function, and since it is antisymmetric, this flow simulates the 
collision of two fluid jets coming from infinite with opposite directions. The 
collision point is the origin, and the (quasi-)shock is generated when two fluid 
particles coUide there. 

We can estabhsh the boundedness of the velocity using the properties of the 
space of functions of bounded mean oscillation (BMO) (for the basic properties 
of this space see, for instance, [24]): 

|k||BMO<C||^||ioc, (20) 

for a finite positive constant C, and 

II^'IIbmo = sup-^ / \v-VQ\dx, (21) 
Q \Q\ JQ 

where Q is any closed interval of M, and vq = Jq vdx. From the very definition 
of the BMO norm we see that for any e > 0, the inequahty 

ll-yllBMO > ;r / \v-v^\dx (22) 

holds, and now we can use the fact that v is an odd function, and thus — 
jt^ vdx — 0, leading us to claim that 

||^^||BMO>^y \v\dx. (23) 

This imphes the estimate 

^ r \v\dx < CII^IIloc, (24) 
2e J—e 
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homogeneous in e, so we can take the hmit e — > to get 

|^;(0,i)| < CII^IIloo, (25) 
where we have used the fact that [25] 

In conclusion, we have shown that while the velocity remains bounded at 
the origin, its first derivative goes to infinity in finite time, or what is the 
same, the fiuid develops either a shock-wave or quasi-shock-wave. Finally, we 
can show that these singular solutions have finite kinetic energy. We know 
that ||^||l2(^) < ||^||l2(0), as proven in Ref.[22]. Using the Calderon-Zygmund 
inequality [26] we can claim that 

\\v\\L2{t) = \\H{9)\\L2{t) < C\\9\\L.{t) < 6-11^11^,2(0), (27) 

for a finite positive constant C, and thus we see that the kinetic energy remains 
bounded for all times. 

To finish our analysis of this model let us point out one very interesting feature 
of it. The particle trajectories are defined, in this case, by 

'i=^=-m. (28) 

what implies, for a particle starting at the origin, that 

for all times, due to the antisymmetric character of H{9). So one particle 
initially located at the origin will stay there for all times. Since this particle 
is stopped, not only the velocity but also the acceleration will vanish for all 
times, implying 

— ^dtv + vd^v = 0. (30) 
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This is, the velocity obeys Burgers equation locally at the origin. This fact is 
very important, since it shows that the pressure vanishes at the origin, and the 
shock is due to an artificial noninteracting character of the particles. These 
considerations push us to conjecture that the shock-wave developed in this 
type of flow is exactly the same as that formed in the Burgers equation. 

4 The hypoviscous Burgers equation 

While Eq.(19) with /x = develops (quasi-)shock- waves in finite time, it has 
been proven that adding a hypoviscosity to this equation results in global 
existence in time of the solution [22]. The corresponding hypoviscous equation 
reads [22] 



where ^ > is the viscosity, 1 < a < 2, and A" is a fractional derivative of 
the Riesz type that is defined from its Fourier transform 



To confirm the analogue between this model and the Burgers equation we 
will analyze the hypoviscous version of the second one, and we will prove 
global existence in time of the solutions. The effect of hypoviscosity on Burgers 
dynamics has been already studied in Ref. [27]; however, in this reference it is 
analyzed the effect of the hypoviscous dissipativity on the Burgers Markovian 
Random Coupling Model, and we will analyze its effect directly on the Burgers 
equation. 

The hypoviscous Burgers equation is 



dte - H{e)dJ = M"^, 



(31) 



(A"^;= -\k\'^e. 



(32) 



dtv — —vVx + jJ.h.'^v. 



(33) 
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Obviously, the case a — 2 corresponds to the standard viscous Burgers equa- 
tion. In this case we will prove the nonexistence of the shock-wave or quasi- 
shock- wave for a situation similar to that of Eq.(19), this is, for two colhding 
fluid jets. We will suppose therefore that t'a; < to simulate this scenario, 
and as assymptotic conditions we will assume that is constant as |x| ^ oo. 
In the following proof, the variables named as C, C", and C denote arbitrary 
constants and its value may change from line to line. 

Deriving once Eq.(33) with respect to x we find 

dtVx = -{vv:,):, + l^A-^v^. (34) 

Now we want to calculate the norm of the velocity 

\vx\\l^ = - j dtv^dx = j {vvx)xdx + fi j A'^Vxdx ^ 0, (35) 



d 

~dr 



where we have used our basic assumption < 0. So we have obtained the 
conservation in time of the norm of Vx- With respect to the norm we 
have 



2di^ 



= - j Vx{vvx)xdx - fi\\A"/\x\\h. (36) 



Leibniz's rule, integration by parts and the application of the boundary con- 
ditions yield the two following equalities 

- j Vx{vvx)xdx = j VxxVxvdx, (37) 



- j Vx{vvx)xdx = - j vldx - j VxxVxvdx, (38) 
that can be combined to provide us with a new formulation of Eq.(36) 

^|lklli. = ^lklli3-/^ir/Vlli.. (39) 
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The third moment of Vx might be estimated as follows 



I I 1 3 II I 1 2 II II 

I'yxllLS < I I'yxl |loo| Pxl 



(40) 



now choose x ^ (1) ct), and use a Sobolev embedding to find 



I LOO 



(41) 



We can now use the Fourier transform of Vr 



(42) 



to claim that 



l|A^/V||i. = ll(A^/V)1li. = J m{vx)f dk^ 
\m(vMk+ f -]^\(vMk< 

k\<R^ ' ""^ J\k\>R 



where we have used the isometry of the Fourier transform in . We still need 



to estimate the second moment of Vx'- 



-\\v. 



\\Vx\\h < Ika^llLllkxIU- < ^ll^^rrllil + 

<^j\vx\\h + c'- (Ii^.iii. + IIA'^/Vlli.), 



(44) 



where we have used a Sobolev embedding. Selecting e small enough we are led 
to conclude 



vx\\h<{i-c'-y' (^iHii.+c'|iiA"/viii.). 



27 V2e 

This inequality, in addition to Eq.(36) yields 



(45) 



X 



-IMh + c-WA'^/^xWh) + 

2e" 2" J 



lA^/^lli. 



(46) 



14 



and employing Eq.(43), choosing a sufficiently large it! and a sufficiently small 
e, we arrive at the desired estimate 



\K\\l^<C. (47) 

Applying a second spatial derivative over Eq.(33) we obtain 

dtVxx = -{vvx)xx + M"^xx- (48) 
We can now compute the LP' norm of v^x 

- j Vxx{vVx)xxdx - ljL\\A'^^^Vxx\\h, (49) 

and by reiteratively using Leibniz's rule, integration by parts, and the bound- 
ary conditions we find that this equation reduces to 



\vxx\\h = Vxvl^dx - i^\\A'^^'^Vxx\\l2. (50) 



The first integral in the right hand side of this equation may be estimated as 
follows 



/ '^Ix'^xdx < \ \Vxx\\l'^\\Vxx\\lA\'^x\W < 



II l|2i || l|2||||2 
^W^xxWl^ "t Il^I I'^xI \'^^) 

and we might continue this chain of inequalities by means of the Sobolev 
embedding 

||^;a:.||ioo < CiWVxxWh + ||A"/V.||iO, (52) 

to conclude 

/e 1 
vlxVxdx < -CdlVxxWh + \\A"'^'^Vxx\\l2) + ^\\Vxx\\lAMh- (53) 
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Choosing e small enough and substituting this result in Eq.(33) we find 

^WVxxWh < C Q + ^IKWh'^ \\Vxx\\h: (54) 



which yields 

\\v:cx\\h<Ce^'. (55) 
We can finish using this result in addition to the Sobolev inequality 

||^^a;||ioo < C(||V^||^2 + ||^^xa;||i2), (56) 

to get the desired estimate 

I^xIIl- <C(l + e^*), (57) 

which prohibits the formation of shock-waves in finite time for the hypoviscous 
Burgers equation. 

5 News from two- and three-dimensional incompressible flows 

In order to understand the physical meaning of the possible blow-ups ap- 
pearing in the Euler equations it is useful to take a look to two- and three- 
dimensional incompressible flow. This is very important to delucidate whether 
or not the one- dimensional models have physical meaning. 

The problem of the blow-up for an incompressible fluid was studied in Ref. [28]). 
They studied four different equations: the Quasi-Geostrophic equation, the 
ideal two-dimensional Magneto- Hydrodynamics equation, the two-dimensional 
Euler equation, and the Boussinesq equation. They found that for all these 
equations, two arcs moving with the fluid cannot collapse in flnite time into 
one single arc, provided the velocity is bounded. This is a strong result against 
the possible formation of (quasi-)shock-waves within this type of flow. Even 
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in three-dimensional incompressible ideal flow there is a similar result proven 
in this same direction [29]. 

As we have seen, a (quasi-) shock- wave develops when two fluid particles collide 
in finite time. Two arcs moving with the fluid are level sets of fluid particles, 
its mutual collapse would imply the formation of a (quasi-) shock- wave. The 
nonexistence of this type of singularity for bounded velocities suggests that 
the possible blow-up of the solution to the three-dimensional incompressible 
Euler equations is related to the breakdown of the nonrelativistic description 
of the fluid, while it has no relation with the phenomenon of turbulence. 
Furthermore, this fact indicates that the appearance of (quasi-) shock- waves 
in one-dimensional models is more related to a mathematical artifact than to 
a real physical phenomenon. 

6 Conclusions 

In this work we have shown that one-dimensional models for three-dimensional 
incompressible hydrodynamics showing the appearance of singular behaviour 
fail to reproduce some of the most important features of the flow. In the best 
cases, the appearance of a (quasi-)shock-wave is due to the strong compress- 
ible character of the flow; actually, we arrive at the absurd conclusion that 
these models for incompressible hydrodynamics are more compressible than a 
compressible flow. This is inferred from the fact that the shock develops due to 
the noninteracting character of the fluid particles in its neighbourhood, what 
produces collisions among them. This strongly suggests that one-dimensional 
settings are too limiting to correctly represent incompressible hydrodynamics. 

On the other hand, the results involving two-dimensional incompressible flow 
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show that the appearance of singularities is related with a divergence of the 
velocity, what would rule out the possibility of (quasi-)shock-wave formation 
within the flow. This would show that possible singularities have nothing to 
say in the transition to turbulence in incompressible flows. Furthermore, the 
divergence of the velocity would indicate the breakdown of the nonrelativistic 
description of the fluid, but would say nothing about the small scale properties 
of it. This indicates that a blow-up of the velocity would have relevance in the 
study of some special fluids, such as interstellar plasmas, which reach velocities 
that may be comparable to that of the light. However, the daily observed 
phenomenon of turbulence, which appears in common fluids at velocities not 
comparable with light speed, seems to be not related with the existence of 
singular solutions to the three-dimensional incompressible Euler equations. 
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